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The muffin-tin model of an effective-mass electron interacting with magnetic ions in semicon- 
ductors is extended to incorporate electrostatic potentials that are present in the case of Mn-based 
III-V compounds (Gai-^Mn^N, Gai-^Mn^As). Since the conduction band electron is repelled from 
negatively charged magnetic ions and attracted by compensating donors, the apparent value of the 
s-d exchange coupling Noct is reduced. It is shown that the magnitude of this effect increases when x 
diminishes. Our model may explain an unusual behavior of electron spin splitting observed recently 
in those two materials in the Mn concentration range x < 0.2%. 



Owing to the possibility of a gradual incorporation of 
magnetism to the well-known semiconductor matrices, 
diluted magnetic semiconductors (DMS) 1 ^^ offer un- 
precedented opportunity for examining energies charac- 
terizing spin dependent couplings between the band car- 
riers and electrons localized in the open magnetic shells. 
Surprisingly, however, a series of recent experiments on 
(III,Mn)V DMS points to our limited understanding of 
the s-d exchange interaction in this important material 
family. The determined s-d exchange integral appears 
to have much smaller magnitude and even opposite 
sign 6 - to that expected according to the present knowl- 
edge on the origin of the s-d coupling in tetrahedrally 
coordinated DMS. 

In this paper, we list first a number of obstacles mak- 
ing a quantitative determination of the exchange inte- 
grals in III-V DMS difficult. We then analyze an addi- 
tional ingredient of these systems, namely the presence 
of Coulomb potentials centered on the magnetic ions as 
well as on compensating donors. We evaluate electron 
wave function in the field of negatively charge magnetic 
impurities and show that the Coulomb repulsion reduces 
the apparent magnitude of the s-d exchange integral. Im- 
portantly, the effect increases with lowering magnetic ion 
concentration x, and becomes particularly significant in 
the experimentally relevant range, x < 0.2%. 4,5,6 

In the case of archetypical II- VI DMS such as 
(Cd,Mn)Te, the ferromagnetic exchange interaction be- 
tween the conduction band electrons and Mn spins is 
described by Noa « 0.2 eV, where No is the cation con- 
centration and a is the s-d exchange integral. This value 
of Noa is about two times smaller than that describing 
the ferromagnetic exchange interaction between the 4s 
and 3d electrons in the free Mn +1 ioni This reduction is 
caused by matrix polarizability and the fact that not only 
cation but also anion s-type wave functions contribute to 
the Bloch amplitude of the conduction band electrons. In 
the case of the valence band holes, the exchange energy 
results from the symmetry- allowed p-d hybridization, the 



typical value of the exchange energy being \Nq0\ ~ 1 eV. 
Within the molecular-field (MFA) and virtual crystal ap- 
proximations (VCA), the exchange spin-splitting of the 
two-fold degenerate conduction and four-fold degener- 
ate valence band is then, s z aM/gfiB and j z (3M/g/j,B, 
where s z = ±1/2 and j z = ±1/2; ±3/2, respectively, 
M = M(T,H) is spin magnetization of the substitu- 
tional magnetic ions characterized by the Lande factor g. 
The proportionality between exchange splittings and in- 
dependently measured magnetization has been demon- 
strated by a variety of magnetooptical and magnetotrans- 
port experiments, and has made it possible to determine 
accurately the values of Noa and No (3 for a number of 
systems 

However, the above simple scenario has been called 
into question in several important cases. First, the or- 
bital and carrier contribution to the measured M has 
to be taken into account^ Second, when the exchange 
energy \Nq[3\ becomes comparable to the valence band 
width, the MFA and VCA break down, particularly in the 
range of small magnetic ion concentrations^S*ii Third, 
the magnitude and sign of (3 depend on the relative po- 
sition of the p and d statesAiS If, therefore, the charge 
state and thus the energy of the relevant d levels can be 
altered by the position of the Fermi energy, the charac- 
ter of p-d exchange will cease to be universal in a given 
material but instead will depend on the doping type and 
magnitude. Fourth, the intensity, and even the sign of 
the magnetic circular dichroism is strongly affected by 
the Moss-Burstein effect. Accordingly, a simple relation 
between positions of the absorption edge for two circular 
light polarizations and the splitting of the bands breaks 
down in the presence of the delocalized or weakly lo- 
calized carrier liquidAi^ This may account for the sign 
reversal of the apparent f3 on going from n-type to p- 
type (Ga,Mn)As>i£ Finally, spin-orbit interactions and 
k-p mixing between bands make spin-splitting away from 
band extrema to be a complex non-linear function of aM 
and /3M as well as of relevant k ■ p parameters. This, in 
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particular, has precluded a conclusive determination of 
the values of the sp-d exchange integrals for narrow-gap 
DMS of mercury and lead chalcogenidesi^i Such multi- 
band effects are especially important in quantum struc- 
tures, where dimensional quantization enhances the ki- 
netic energy of the carriers and, thus, the effects of the 
k ■ p coupling. Indeed, an anomalous behavior of elec- 
tron spin-splitting in quantum wells of (Cd,Mn)Te and 
(Ga,Mn)As has been assigned to the k ■ p admixture of 
the valence band states to the electron wave function&i^ 

Despite the difficulties in the precise determination of 
the exchange integrals, particularly in quantum struc- 
tures and systems containing carriers, a series of recent 
experiments suggesting anomalous magnitude and sign of 
a in (III,Mn)V DMSte& call for a detail consideration. 
In particular, Heimbrodt et alX detected spin-flip Raman 
scattering of conduction band electrons in Gai-^Mn^As, 
and evaluated Nqcx — 23 meV for x = 0.1%. Even a 
lower value \Noa = 14 ± 4| meV was found by Wolos 
et al.^ who analyzed the broadening by the electrons 
of the Mn spin resonance line in n-Gai-zMn^N with 
0.01% < x < 0.2%. More recently, Myers et al& ex- 
amined spin precession of the electrons in Gai-^Mn^As 
quantum wells of the thickness between 3 and 10 nm, 
and Mn content x up to 0.03%. As a result of afore- 
mentioned admixture of the valence band states, the ob- 
served sign of the exchange splitting is negative. The 
value Nqci = — 90 ± 30 meV was determined under a sim- 
plified assumption that the spin-splitting is proportional 
to magnetization, and by extrapolating the resulting ap- 
parent exchange energy NqQ! to the infinite quantum well 
width.- In contrast to the striking finding listed above, 
a large positive value Nga « 0.5 eV is consistent with 
intraband magnetoabsorption in Ti-Ini-^Mn^As with a 
relatively high Mn content, x > 2.5%^ 

We point here to an additional mechanism that may 
contribute to the anomalous behavior of electron-spin 
splitting in III-V DMS. We note that the electric charge 
of the Mn 2+ ion replacing e.g. a Ga 3+ ion in the lattice 
of a III-V compound (like GaN and GaAs) is a source 
of a repulsive electrostatic potential. Furthermore, the 
studied samples are either n-type^ or at least highly com- 
pensated, as evidenced by the presence of electron spin- 
flip Raman scattering^ and donor-related luminescence ^ 
This indicates the existence of attractive potentials as- 
sociated with ionized non-magnetic donors. Thus, the 
probability of finding a conduction band electron at the 
core of the magnetic ion is reduced, and hence the ap- 
parent value of the exchange energy (the observed spin 
splitting) is diminished. It worth noting that a possi- 
bility that the Coulomb potentials could affect the ap- 
parent value of the exchange integrals has already been 
mentioned in the context of divalent Mn in GaN^ and 
trivalent Fe in rlgSe^ 

To evaluate a lower limit of the effect we neglect the 
presence of compensating donors and calculate the ap- 
parent s-d exchange integral a ap for an electron subject 
to the repulsive potential generated by the Mn acceptors. 
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FIG. 1: Wave functions of spin-up and spin-down carriers (the 
Coulomb term included) for Gai-^Mn^N and R = 75a.u. 
(x = 0.0087%). 
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FIG. 2: Wave functions of spin-up and spin-down carriers (the 
Coulomb term included) for G&i- x }An x As and R — 250 a.u. 
(x = 0.00047%). 



We follow a Wigner-Seitz-type approach put forward by 
one of us and co- workers^ to describe the interaction of 
the carrier spin with the Mn ions in the case of the strong 
coupling limit, that is when the depth of the local Mn po- 
tential is comparable to the carrier band width. It has 
been found in the subsequent worksii that the correc- 
tions to the Wigner-Seitz approach caused by a random 
distribution of Mn ions are quantitatively unimportant. 

We consider a Mn ion with the 5/2 spin Si located 
at which interacts with the carrier via the Heisen- 
berg term I(r — Ri)s ■ Si. The form of the func- 
tion I(f — Ri) makes the interaction local: it van- 
ishes outside the core of the Mn ion. For simplicity, 
I(f — Ri) = a9{b — \ f — Ri\). The exchange energy is 
then a = f d 3 fl(r) — a ■ §7r& 3 - Moreover, in case of 
III-V compounds considered here, the impurity gener- 
ates an electrostatic potential. If screening by the elec- 
trons is present, as in case of n-Gai^^Mn^N, this po- 
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FIG. 3: The assumed dependence e = e(r). 



tential is e 2 exp(— Xr) / (iweosr) , where e is the static di- 
electric constant, and the screening parameter A is given 
by A 2 = e 2 Af(£ F )/(e e), where M{£ F ) = ^n/kT F il For 



the Gai_ x Mn x N samples," 



10 cm 3 corresponds to 



T F w 890 K (£ F w 0.12 eV), and therefore 1/A « 1.6 nm. 

In the spirit of the Wigner-Seitz approach we assume 
that the carrier energy E and the envelope function ip(r) 
are given by the ground state s solution of the one-band 
effective mass equation which contains the potential U(r) 
created by the magnetic ion located at r = 0. The 
standard one-impurity boundary condition ij)(r) — ► for 
r — ► oo is replaced by the matching condition ip'(r) = 
at r = R to take into account the presence of other mag- 
netic ions. The value R is determined by the concen- 
tration of the magnetic ions x according to the equation 
(47ri? 3 /3) _1 = Nqx. The exchange interaction is mod- 
elled by a square well potential U6(b — r) superimposed 
on the electrostatic potential of an elementary charge lo- 
cated at r = 0. The potential U = ±|a is, of course, 
different for spin-down and spin- up carriers. 

We first ignore free carrier screening, A — > 0. The 
solution of the time-independent Schrodinger equation 
for the conduction band electron is then 

i>(r) = c Q exp(-/3r)$(l + -; 2; 2/3r) = c a f (1) 

for < r < b, and the following linear combination for 
b <r < R 



V>(r) = ciexp(-/3'r)*(l + ^;2;2/3V) + 

+ c 2 exp(/3V)*(l-^;2;-2/?V) 
= cig + c 2 h, 



(2) 



where A = e 2 m* /(Aneoeh 2 ), (3 = [2m* (U- E)]? /%, /3' = 
[2m*(-£)]5/S (notice that changing the sign of (3 leaves 
tjj invariant, while changing the sign of (3' interchanges 
ci with ci\ also, $ and 4" are not in general linearly 
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FIG. 4: The dependence of the ratio of the apparent and bare 
exchange energies a on x for for Gai_a;Mna;N and various 
models of screening. 



independent). We used the symbols $, 9 for the con- 
fluent hypergeometric functions ii*i(a;6; z), U(a;b; z)w& 
The constants cq, Ci, ci are determined by the continu- 
ity conditions ip(b~) = ?p(b + ), ip'(b~) = ip'(b + ). Solving 
those two equations we obtain an equation for E, 



w f , h {b)g'{R)-w f , g (b)h'{R) 

Wg,h{b) 



= 0, 



(3) 



In 



where by Wf t3 we denoted the Wronskian fg' — f'g 
the following, ip(r) is normalized as ip(0) = c = 1. 

We assume the following parameters for Gai-^Mn^N: 
m* = 0.22 m e , N = 4.38- 10 22 cm" 3 = 0.006495 a.u., e = 
8.9; and the following for Gai-^Mn^As: to* = 0.067 m e , 
N = 2.21 • 10 22 cm- 3 = 0.003281 a.u., e = 12.9. In the 
experiments, samples were used with 0.01% < x < 0.2% 
of Mn in GaN£ and with 0.0006% < x < 0.03% of Mn 
in GaAs£ Those concentrations correspond to R up to 
about 75 a.u. for GaN and up to about 250 a.u. for GaAs. 

To visualize the effect of the Coulomb term in the 
Mn potential, we have calculated the energies and 
wave functions including and disregarding the additional 
Coulomb term for both GaN (b — 2 a.u. « 0.1 nm, 
a = 0.0371 a.u. = 1.0 eV) and GaAs (6 = 2 a.u. w 0.1 nm, 
a = 0.0735 a.u. = 2.0 eV). These parameters correspond 
to N a = 0.22 eV, a value for CdSi& We have found that 
when calculating a ap /a, the details of the exchange po- 
tential (like the values of b and a within the expected 
range) are not quantitatively important. 

In order to take into account the fact that the core and 
lattice polarizability decrease at small distances, e — > 1 
for r — > 0, we interpolate e(r) between e(0) = 1 and 
the macroscopic value attained at a distance of the bond 
length. The assumed dependence, presented in Fig. |3 is 
similar to that of the Thomas-Fermi models When e = 
e(r) and/or free carrier screening is included, we find the 
solution ijj(r) of the Schrodinger equation for the given 
potential U(r) numerically, as Eqs. Q and @ are only 
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valid for the Coulomb potential. Then, the spin splitting 
for a given value of x (or for the corresponding R) is 
evaluated as the difference of the energy E calculated 
for the spin-up and spin-down carriers from the equation 
i/j'(R) = 0. Here, ip(r) is the numerical solution of the 
Schrodinger equation with the potential that is different 
for spin- up and spin-down carriers. 

The results of our calculations of a ap /a as a function 
of the Mn ion concentration x are presented in Fig. 0] 
(Gai_ x Mn x N) and in Fig. (Gai_ x Mn x As). Indepen- 
dently of assumptions concerning screening, in both ma- 



terials a ap /a diminishes significantly when x decreases, 
up to factor of three in the experimentally relevant range 
of x. However, this reduction of a ap /a is still smaller 
than that seen experimentally^ 5 - presumably because 
of an additional effect coming from the presence of at- 
tractive potentials brought about by compensating non- 
magnetic donors. 

In summary, we have enlisted a number of effects that 
renders an accurate experimental determination of the 
sp-d exchange integrals difficult, particularly in cases 
when both p-like and s-like states contribute to the 
carrier wave function. The interaction of conduction 
band electrons with the magnetic ions in (Gai_ x Mn x N, 
Gai-^Mn^As) has been considered quantitatively taking 
into account the electrostatic potential created by the 
magnetic ion. A substantial reduction in the magnitude 
of the apparent exchange energy has been found at low 
Mn concentrations, and interpreted as coming from the 
decrease of the carrier probability density at the core of 
the magnetic ion caused by the electrostatic repulsion. It 
has been suggested that this effect, enhanced by an at- 
tractive potential of compensating donors, accounts for 
abnormally small values of the exchange spin splitting ob- 
served experimentally in III-V DMS containing a minute 
amount of Mn. 4 ' 5 ' 6 In view of our findings, the presence 
of electrostatic potentials associated with magnetic ions 
makes that the apparent exchange energies should not be 
viewed as universal but rather dependent on the content 
of the magnetic constituent and compensating donors. 
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